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ABSTRACT 

In  April  2005,  a  laser  propagation  experiment  was  conducted  over  a  470m  horizontal  maritime  path.  Scintillation 
measurements  of  a  divergent  Gaussian  beam  wave  were  taken  simultaneously  for  different  receiver  aperture  sizes. 
Terrestrial  scintillation  theory  combined  with  a  numerical  algorithm  was  used  to  infer  the  atmospheric  parameters  Cn2 
and  la  from  the  optical  maritime  scintillation  measurements.  This  paper  presents  the  initial  results. 

1.  INTRODUCTION 

In  April  2005,  a  laser  propagation  experiment  was  conducted  over  a  470m  horizontal  maritime  path,  located  at  the 
Naval  Research  Laboratory’s  (NRL’s)  facilities  at  the  Chesapeake  Bay,  Maryland.  Scintillation  measurements  of  a 
divergent  Gaussian  beam  wave  were  taken  simultaneously  for  different  receiver  aperture  sizes.  Scintillation  theory  for 
Gaussian  beam  waves,  recently  developed  by  Andrews  et  al  for  terrestrial  conditions,1’2  combined  with  a  numerical 
algorithm  was  used  to  infer  the  refractive  index  structure  constant  C„  and  the  inner  scale  of  turbulence  10  from  the 
optical  scintillation  measurements.  In  this  experiment,  the  scintillation  data  was  collected  under  the  conditions  of  weak 
optical  turbulence.  The  proposed  method  of  inferring  atmospheric  parameters  can  be  extended  into  the  moderate-to- 
strong  fluctuation  regime,  in  which  case  values  of  the  outer  scale  of  turbulence  Ln,  as  well  as  C„2  and  ln,  can  be 
determined.3 

This  work  contains  the  initial  results  of  applying  scintillation  theory  developed  for  terrestrial  conditions  to  maritime 
scintillation  data.  The  ultimate  goal  is  to  make  corrections  to  the  scintillation  theory  for  maritime  conditions  and  redo 
this  analysis  to  determine  the  impact  of  such  corrections. 

2.  INFERRING  ATMOSPHERIC  PARAMETERS  FROM  OPTICAL  MEASUREMENTS 

The  goal  of  this  work  is  to  infer  the  refractive  index  structure  parameter,  C„2,  and  the  inner  scale  of  turbulence,  /#,  from 
optical  scintillation  measurements  taken  under  weak  turbulence  conditions  in  a  marine  environment.  The  scintillation 
index,  defined  as  the  normalized  variance  of  the  irradiance, 
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where  /  represents  the  irradiance  of  the  optical  wave  and  <•>  denotes  ensemble  average,  is  known  to  be  sensitive  to  both 
C„ 2  and  /o  under  all  conditions  of  optical  turbulence.  The  outer  scale  of  turbulence,  L0,  has  a  negligible  effect  on  the 
scintillation  index  under  weak  irradiance  fluctuations,  but  its  influence  begins  to  emerge  as  the  strength  of  the  irradiance 
fluctuations  increases.  Hence,  scintillation  measurements  taken  in  the  weak  optical  turbulence  regime  cannot  be  used  to 
infer  the  outer  scale. 
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When  comparing  theory  with  experimental  data  obtained  with  larger  apertures,  aperture  averaging  effects  (a  reduction 
in  scintillation  with  increasing  receiver  aperture  size)  have  to  be  taken  into  account.  As  a  result,  by  ignoring  outer  scale 
effects,  the  theoretical  expressions  for  the  scintillation  index,  at  a  fixed  range  and  wavelength,  depend  on  three 
parameters,  C„2,  l0  and  D, 


°7  (A)>  ,D)  -  exp  crlnv(/0,C„,Z))  +  cjln>,(/{l,C„,Z))  1, 
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where  D  denotes  the  diameter  of  the  receiver  aperture,  and  (7lnv  and  CTln  are  the  large  scale  and  small  scale  log- 

irradiance  variances,  respectively.  Expressions  for  the  log-irradiance  variances  are  provided  in  Appendix  A.  By 
specifying  two  different  collecting  aperture  sizes,  D,  and  D2,  and  measuring  the  corresponding  scintillation  index 
experimentally,  <r;2exp ,  the  result  is  a  system  of  two  nonlinear  equations  with  two  unknowns, 

cjj  (l0 ,  C2„ ,  Dj )  =  alexp  ( Dj ),  j  =  1, 2  ,  (3) 

where  l0  and  Cn2  are  the  unknowns  to  be  determined  and  CTj  ,  D )  is  the  theoretical  scintillation  index  defined  in 

(2). 

Due  to  the  complexity  and  nonlinearity  of  (3),  the  system  was  solved  numerically.  However,  rather  than  solving  the 
system  of  equations  in  (3)  directly,  the  problem  was  rewritten  as  a  minimization  problem  by  introducing  the  function 

m , C2n )  =  t \aj  (/„ , C2n  ,Dj)~  (Ty  exp (Dj  )T  .  (4) 

j= i 

A  customized  Mathematica  program  based  on  the  downhill  Simplex  method4'5  was  used  to  find  the  I0  and  Cn2  values 
that  minimize  f[lo,C„2)  in  (4).  The  downhill  Simplex  method  is  a  multidimensional  minimization  algorithm  that  requires 
only  function  evaluations,  not  derivatives.  The  simplex  is  a  geometrical  figure  consisting,  in  N  dimensions,  of  N+l 
points  (vertices)  and  their  interconnecting  line  segments.  In  two  dimensions,  the  simplex  is  a  triangle.4 

In  addition  to  an  initial  estimate,  the  downhill  Simplex  method  requires  N  extra  points  to  form  the  initial  simplex.  These 
points  are  determined  by  specifying  the  initial  estimate,  Po,  and  an  initial  step-size,  8.  Both  P0  and  8  are  N-dimensional 
vectors.  The  additional  N  points  are  then  given  by  P;  =P0+Sj  e, ,  where  8;  is  the  ith  element  in  the  step-size  vector  8  and  e; 
is  a  unit  vector.  The  algorithm  moves  the  simplex  downhill  on  the  N-dimensional  surface  until  it  reaches  a  minimum, 
where  the  simplex  contracts  itself  in  all  directions.  The  lowest  point  of  the  simplex  is  chosen  as  the  solution  of  the 
minimization  problem.4 

The  termination  criterion  for  the  Simplex  method  cannot  be  chosen  to  require  a  certain  tolerance  for  a  single 
independent  variable.  The  algorithm  can  be  programmed  to  terminate  when  either  the  decrease  in  the  function  value  in 
one  step  is  smaller  than  some  tolerance,  or  the  vector  distance  moved  in  one  step  becomes  smaller  than  some  tolerance. 
Due  to  the  experimental  error,  the  minimum  value  of  (4)  will  most  likely  never  be  exactly  zero.  Therefore,  it  is  better  to 
use  a  termination  criterion  that  depends  on  the  variables,  rather  than  the  function  value.  Since  the  simplex  contracts 
itself  around  a  minimum,  the  “simplex  size”,  defined  as  the  sum  of  the  Euclidean  vector  distance  between  all  vertices  in 
the  simplex,  was  used  as  a  termination  criterion. 


3.  EXPERIMENT 


3.1  Experimental  setup 

The  experiment  was  conducted  on  April  6,  2005,  at  NRL’s  facilities  at  the  Chesapeake  Bay,  Maryland.  A  continuous 
wave  laser  operating  at  1550nm  was  propagated  horizontally  470m  over  open  water  between  two  piers.  The  beam 
divergence  was  estimated  to  be  3  mrad  (full  angle).  The  beam  was  detected  by  three  germanium  photodiodes  of  1mm, 
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5mm  and  13mm  diameters  (Judson  J16TE2-8A6-R01M,  J16TE2-8A6-R05M,  J16TE1-P6-R13M-HS),  respectively, 
and  a  12.7cm  diameter  telescope  (Maksutov-Cassegrain)  combined  with  a  germanium  detector.  To  reduce  background 
radiation,  two  1550nm  notch  fdters  were  placed  in  front  of  the  5mm  and  13mm  Ge  detectors,  and  their  field  of  view  was 
limited  by  tubes  of  2.54cm  diameter  and  45.72cm  length.  In  order  to  receive  signal  on  the  1mm  Ge  detector,  no  filter 
was  used  and  the  length  of  tube  was  reduced  to  15.24cm.  The  telescope  detector  had  neutral  density  filters  placed  in 
front  of  it. 


Two  tripods  were  used  to  hold  the  detectors,  one  with  the  3  Ge-detectors  and  one  with  the  telescope.  All  detectors  were 
positioned  about  3 -4m  above  the  water  surface.  Each  of  the  three  Ge-detectors  was  connected  to  a  pre-amplifier  (Judson 
PA-7-70),  and  the  voltage  signal  was  captured  with  a  sampling  rate  of  10kHz  by  a  data  acquisition  card  (National 
Instillments  PCI-607  IE)  and  a  custom-written  Lab  VIEW  software.  Data  was  taken  in  5  minute  intervals  and  the 
background  radiation  was  measured  either  before  or  after  each  run.  During  the  data  collection,  the  weather  was  sunny 
and  very  windy,  with  air  temperatures  in  the  range  15-20  °C. 

3.2  Data  analysis 

A  Matlab  program  was  written  to  calculate  the  mean  signal  and  the  scintillation  index  for  each  detector  and  each  minute 
of  collected  data.  In  order  to  obtain  the  true  experimental  scintillation  index  (SI),  cr/exp,  the  background  noise  was 
subtracted  from  each  data  point.  The  background  noise,  VBG  ,  was  estimated  to  be  the  mean  signal  of  the  background 
measurement  taken  either  before  or  after  each  5-minute  interval  of  collected  data.  Rather  than  subtracting  VBG  from 
each  data  point  and  recalculating  the  SI  from  the  corrected  raw  data,  the  corrected  SI  can  be  calculated  according  to 


a 


2 

I, ex  p 


(5) 


where  ),c/data  and  Fdata  are  the  SI  and  mean  signal  (voltage)  calculated  by  the 
Matlab  program  for  each  minute  of  collected  data.  The  average  signal-to-noise  ratio 
(SNR)  for  all  data  runs  varied  for  the  different  detectors  and  is  shown  in  Table  1.  For  the 
5mm,  13mm  and  12.7  cm  apertures,  the  background  correction  resulted  in  a  minor 
change  in  the  scintillation  index,  making  cr/2exp  on  average  only  0.2%,  2%  and  0.06% 
higher  than  of  data  ,  respectively.  However,  for  the  1mm  aperture,  cr/2exp  was  more  than 
50%  higher  than  cr/datll.  In  addition,  cr/exp  for  the  1mm  aperture  was  consistently  lower 
than  cr/exp  for  the  5mm  aperture,  which  is  contradicting  aperture  averaging  theory. 

Hence,  the  data  from  the  1mm  could  not  be  used. 

Before  applying  the  downhill  Simplex  method  to  minimize  ftJo,Cn2)  in  (4),  the  variables 
l0  and  C„:  had  to  be  made  comparable  in  size.  Hence,  I0  was  specified  in  mm  and  C„2  in 
m'2,J  on  a  logi0-scale.  A  general  two  dimensional  vector  was  set  to  V  =  ( l0,  logi0[C„2  ]  ), 
the  initial  step-size  was  set  to  8  =  (  1,  0.1  )  and  the  initial  estimate  was  chosen  as  P0  =  (  3,  log10[C„2jnjtiai  ]  ).  V,  P0  and  5 
all  have  units  (  mm,  logio[m"23]  ).  C„2jnj,jai  was  determined  from  the  relation 

=  l+CTtap)D7/3L'3  (6) 

where  D  is  the  receiver  aperture  diameter,  L  the  propagation  distance  and  ),0/2exp  is  the  experimental  scintillation  index 
associated  with  D.  Eq.  (6)  provides  a  good  estimate  of  C„  for  large  apertures,  so  the  data  from  the  12.7cm  aperture  was 
used  to  calculate  C„2jnjtja|.  The  termination  criterion  was  chosen  as  the  “simplex  size”  <  10"9.  This  value  was  chosen  to 
ensure  a  unique  solution. 

A  customized  Mathematica  program  employing  the  downhill  Simplex  method  was  used  to  perform  a  minimization  for 
each  minute  of  data.  Since  scintillation  data  was  available  from  three  different  receiver  apertures  (5mm,  13mm  and 


Table  1. 


Aperture 

diameter 

SNR  [dB] 

1  mm 

5.9 

5  mm 

30.5 

13  mm 

20.0 

12.7  cm 

35.6 
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12.7cm  diameter),  but  only  two  are  required  for  the  minimization,  three  different  combinations  of  aperture  sizes  could 
be  used  to  infer  the  atmospheric  parameters.  The  numerical  scheme  was  tested  with  simulation  data,  which  showed  that 
the  solutions  obtained  from  the  5mm  and  13mm  apertures  were  unstable,  i.e.  very  sensitive  to  small  changes  in  the 
scintillation  data.  Hence,  this  combination  of  aperture  sizes  was  disregarded.  Results  for  C„  and  inner  scale  were 
obtained  with  both  the  5mm  and  12.7cm  apertures  and  the  13mm  and  12.7cm  apertures.  In  both  cases,  the  separation  in 
aperture  size  was  large  enough  to  ensure  stable  solutions  of  C„~  and  l0. 

4.  RESULTS  AND  DISCUSSIONS 

The  downhill  Simplex  method  converges  to  a  unique  solution  of  the  atmospheric  parameters  ( C„2  and  In)  which  agree 
with  our  physical  intuition,  i.e.  a  C„2-value  in  the  range  10'14  -  10'13  m" 2/3  and  In  of  the  size  of  a  few  millimeters.  The 
inferred  Cn2-  and  10- values  are  displayed  in  Figs.  1  and  2,  respectively,  as  a  function  of  time  of  day.  The  displayed  data 
is  from  April  6,  2005  and  the  results  are  shown  for  both  the  13mm  and  12.7cm  apertures  (gray)  and  the  5mm  and 
12.7cm  apertures  (black).  The  majority  of  the  inferred  /^-values  lie  between  1.5mm  and  2.5mm.  These  low  values  of  the 
inner  scale  of  turbulence  can  be  explained  by  the  extremely  windy  conditions  during  the  experiment,  since  the  inner 
scale  is  known  to  be  inversely  related  to  the  average  wind  speed.  For  some  minutes  of  data,  the  numerical  scheme 
resulted  in  inner  scale  values  very  close  to  zero  and  very  high  C„2-  results  (on  the  order  of  10"l:  m" 2/3).  These  results  do 
not  make  physical  sense  and  should  be  disregarded.  (The  high  Cn2-  values  are  not  shown  in  Fig.  1  since  they  lie  outside 
the  range  of  the  majority  of  the  data). 

Fig.  3  shows  a  comparison  of  the  inferred  atmospheric  parameters  obtained  with  different  combinations  of  receiver 
aperture  sizes  (13mm  and  12.7cm  apertures  and  the  5mm  and  12.7cm  apertures)  for  each  minute  of  data.  The  few 
minutes  with  an  inner  scale  solution  close  to  zero  have  been  omitted  from  Fig.  3.  As  shown  in  Fig.  3a,  the  C„ 2-  values 
inferred  from  the  5mm  and  12.7cm  apertures  are  consistently  larger  than  those  inferred  from  the  13mm  and  12.7cm 
apertures.  The  absolute  difference  between  the  two  C„2-  values  is  on  average,  for  all  data  runs,  9T0"15  m"23.  As  shown 
in  Fig.  3b,  the  two  different  combinations  of  aperture  sizes  yield  almost  identical  values  of  the  inner  scale  of  turbulence. 
On  average,  for  all  data  runs,  the  absolute  difference  of  the  two  inferred  inner  scale  solutions  is  0.06mm. 

The  accuracy  of  the  solution  to  the  Downhill  Simplex  method  is  related  to  the  success  of  the  minimization  of  f[I0,C„2)  in 
(4),  i.e.  how  small  the  difference  between  the  measured  and  the  theoretical  scintillation  index  is  for  each  aperture  size. 
The  success  of  the  minimization  can  be  measured  in  terms  of  the  error  of  the  theoretical  scintillation  index  relative  to 
the  measured  scintillation  index,  which  is  shown  for  each  aperture  size  in  Fig.  4.  The  theoretical  scintillation  index  is 
calculated  by  inserting  the  inferred  C„2-  and  /«- values  into  (2)  for  each  receiver  aperture  size.  Fig.  4a  shows  the  relative 
error  calculated  with  the  atmospheric  parameters  inferred  from  the  measured  scintillation  data  from  the  5mm  and 
12.7cm  apertures,  while  Fig.  4b  displays  the  same  quantity  obtained  from  the  13mm  and  12.7cm  apertures. 


Fig.  1.  The  inferred  C,2-  values  as  a  function  of  time  Fig.  2.  The  inferred  lo  -  values  as  a  function  of  time  of 
of  day.  day. 
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(b) 


Fig.  3.  Comparison  between  the  inferred  C„-  values  (a)  and  /^-values  (b)  obtained  from  the  experimental  data 
from  either  the  13mm  and  12.7cm  apertures  or  the  5mm  and  12.7cm  apertures. 
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♦  5mm  aperture  o  13mm  aperture  o  12.7cm  aperture 
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Fig.  4.  The  relative  error  of  the  theoretical  scintillation  index  (SI),  compared  to  the  experimental  SI  for  all  three 
aperture  sizes.  The  theoretical  SI  for  all  three  aperture  sizes  were  calculated  with  the  l0-  and  C„-  values  inferred 
from:  (a)  the  5mm  and  12.7  cm  apertures,  (b)  the  13mm  and  12.7cm  apertures. 
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As  shown  in  Fig  4,  the  relative  error  for  the  smallest  aperture  used  in  the  minimization,  the  5mm  aperture  in  Fig  4a  and 
the  13mm  aperture  in  Fig  4b,  is  almost  zero,  on  average  0.08%  and  0.16%,  respectively,  while  the  relative  error  for  the 
12.7cm  aperture  is  significantly  larger  in  both  cases,  on  average  49%  and  51%,  respectively.  The  significant  difference 
in  the  relative  error  for  the  two  apertures  used  in  the  minimization  occurs  since  the  measured  scintillation  index  for  the 
smallest  aperture  is  at  least  one  order  of  magnitude  larger  than  that  for  the  12.7cm  aperture.  Therefore,  the  relative  error 
for  the  smallest  aperture  is  favored  in  the  minimization  scheme  and  becomes  almost  zero.  The  large  relative  error  for  the 
12.7cm  aperture  is  a  concern,  however,  since  aperture  averaging  theory  is  not  exact,  a  large  relative  error  for  the  largest 
aperture  is  acceptable.  The  poor  agreement  between  theory  and  measurement  for  the  large  aperture  case  may  be  due  to 
the  fact  that  the  scintillation  theory  used  was  developed  and  validated  for  terrestrial  conditions,  while  this  experimental 
data  was  collected  in  a  maritime  environment.  Further  investigation  is  needed  in  order  to  determine  if  the  scintillation 
theory  needs  to  be  modified  for  the  maritime  case. 

The  third  relative  error  shown  in  Fig  4  displays  the  results  for  the  aperture  that  was  not  used  in  the  minimization 
scheme,  the  13mm  aperture  in  Fig  4a  and  the  5mm  aperture  in  Fig  4b.  This  error  shows  how  well  the  inferred 
atmospheric  parameters  can  be  used  to  predict  the  scintillation  value  at  an  arbitrary  receiver  aperture,  compared  to  its 
measured  value.  On  average,  this  relative  error  is  11.7%  in  Fig  4a  and  11.8%  in  Fig  4b.  In  this  context  this  is  an 
acceptable  error  and  justifies  the  use  of  this  method  to  infer  C„‘  and  the  inner  scale  of  turbulence.  However,  if  the  theory 
is  adjusted  to  properly  describe  a  maritime  environment,  improvements  of  these  results  may  be  expected. 

5.  SUMMARY 

The  purpose  of  the  experiment,  which  was  conducted  in  April  2005,  was  to  obtain  simultaneous  scintillation 
measurements  from  two  receivers  of  different  aperture  size  and  use  the  results  to  determine  the  atmospheric  parameters 
C„ 2  and  /0  along  a  maritime  propagation  path.  The  parameters  were  inferred  from  a  recently  developed  scintillation 
theory  that  is  believed  to  be  valid  under  weak  to  strong  irradiance  fluctuations.  Receiver  aperture  sizes  of  1mm,  5mm, 
13mm  and  12.7cm  were  placed  at  470m  from  a  laser  source  operating  at  1.55pm,  resulting  in  scintillation  data  from  the 
weak  regime  of  optical  turbulence.  These  aperture  sizes  were  not  considered  an  optimal  choice  for  this  experiment  but 
chosen  primarily  on  the  basis  of  availability. 

Given  the  optical  measurements  from  two  of  the  receiver  apertures,  a  numerical  algorithm  was  developed  in  order  to 
extract  the  atmospheric  parameters  from  the  complex  theoretical  expressions.  Due  to  poor  signal-to-noise  ratio,  the  data 
from  the  1mm  receiver  aperture  was  disregarded.  The  C„2-  and  /0-values  inferred  from  the  scintillation  data  from  the 
5mm  and  12.7cm  apertures  and  the  13mm  and  12.7cm  aperture  were  in  good  agreement.  The  numerical  scheme  was  not 
stable  with  respect  to  small  changes  in  the  experimental  data  for  the  5mm  and  13mm  apertures.  Hence,  the  latter 
combination  of  aperture  sizes  could  not  be  used  to  infer  the  atmospheric  parameters. 

In  conclusion,  scintillation  measurements  from  two  receiver  apertures,  with  a  large  separation  in  aperture  size,  can  be 
used  to  infer  path  average  values  of  C„  and  I0  in  the  weak  regime  of  optical  turbulence  in  a  maritime  environment. 
Since  the  scintillation  theory  used  was  developed  and  validated  for  terrestrial  conditions,  further  investigation  is  needed 
in  order  to  determine  if  modifications  are  required  for  the  maritime  case. 

The  main  advantage  of  the  technique  presented  here  is  a  straightfoiward  extension  into  the  moderate-to-strong 
fluctuation  regime  and  long  propagation  paths.  In  the  moderate-to-strong  fluctuation  regime,  scintillation  data  from 
three  different  receiver  aperture  sizes  are  required  and  the  outer  scale  of  turbulence  L0  is  obtained  in  addition  to  C„2  and 
10.  Existing  methods  of  inferring  inner  scale  values  are  limited  to  weak  fluctuations  conditions  and,  hence,  short 
propagation  paths. 

This  experiment  was  repeated  in  August  2005,  when  scintillation  data  was  collected  from  seven  different  apertures. 
Further  data  analysis  is  in  progress.  Future  work  includes  experiments  over  longer  maritime  propagation  paths,  in  the 
moderate-to-strong  optical  turbulence,  where  path  average  values  of  all  three  atmospheric  parameters  (C„2,  10  and  L0) 
can  be  determined.  In  addition,  further  analysis  on  a  maritime  spectral  model  of  the  irradiance  fluctuations  and  its 
impact  on  scintillation  is  necessary.  If  a  modification  of  the  existing  terrestrial  scintillation  theory  is  required,  this 
analysis  should  be  repeated  and  comparisons  made  to  the  present  results. 
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APPENDIX  A 


At  the  transmitter,  a  Gaussian  beam  wave  is  defined  in  terms  of  the  beam  radius  W0  and  the  phase  front  radius  of 
curvature  F0.  These  parameters  are  related  to  the  (hard)  transmitter  aperture  diameter  d  and  the  half-angle  transmitter 
beam  divergence  6div  ,  according  to 


W0  = 


2^2’ 


F<>=- 


fVn 


ydiv 


(A.1) 


A  propagating  beam  is  characterized  by  the  beam  parameters,  which  in  the  transmitting  plane  are  defined  by, ( 


@o  =  l - , 

^0 


A0  - 


2  L 


(A.2) 


where  L  is  the  propagation  distance  and  k  is  the  optical  wave  number.  At  the  receiver,  located  at  propagation  distance  L, 
the  receiver  beam  parameters  are  defined  by, 6 

@n 


0  = 


®0  +  Aq 


0  =  1-0. 


A  =  A° 


®0  +  Aq  ' 


(A3) 


When  outer  scale  effects  are  ignored,  the  theoretical  expression  for  the  scintillation  index  of  a  receiver  with 

2  7 - A  T\  7 


aperture  diameter  D  is  a  function  of  Cn  ,  Iq  and  D , 

a]  ( /0 ,  Cl ,  D)  =  exp  [o£x  (/0  ,C*,D)  +  o&y  (/0 ,  C„2 ,  £>)]  - 1 , 


(A.4) 


2  2 

where  <r^r  and  cr\nv  are  large-scale  and  small-scale  log-irradiance  scintillation,  respectively.  For  mathematical 
convenience,  the  following  dimensionless  quantities  that  depend  on  C„2,  lo  and  D  are  introduced  7 


aa«)  =  ^.  aUC;)  =  l.23CF116  nc(Z))=16i 


kin 


kD1 


(A- 5) 


where  L  is  the  propagation  distance  and  k  is  the  optical  wave  number.  For  the  on-axis  portion  of  a  Gaussian  beam  wave, 

2  *  7 

(7\nx  is  defined  by 
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where 


rjxd(l(),cfi,D) 
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(A.  6) 


(A.  7) 
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(A.8) 


and  0  and  A  are  Gaussian  beam  parameters  in  the  receiver  plane.  The  small-scale  log-irradiance  scintillation,  <J\nv  , 
for  a  beam  wave  (on-axis),  is  defined  by  7 


aly(l0,C2„,D)  = 
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where 
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and 


<Pl  =  tan 
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where  0  and  A  are  Gaussian  beam  parameters  in  the  receiver  plane. 
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